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Abstract 

In this paper we announce the result of asymptotic dynamics of soHtons of 
nonhnear Schrodinger equations with external potentials. To each local minima 
of the potential there is a soliton centered around it. Under some conditions on 
the nonlinearity, the potential and the datum, we prove that the solution can 
be decomposed into two parts: the soliton and the term dissipating to infinity. 

1 Introduction 

Problem. In this paper we study dynamics of solitons in the generalized non- 
linear Schrodinger equation (NLS) in dimension d ^ 2 with an external potential 
Vh ■■ M" ^ M, 

*^ = -A^ + yw^-/(|^|2)v.. (1) 

Here ft. > is a small parameter giving the length scale of the external potential 
in relation to the length scale of the = solitons (see below) , A is the Laplace 
operator and /(s) is a nonlinearity to be specified later. We normalize /(O) = 0. 
Such equations arise in the theory of Bose-Einstein condensation ^, nonlinear 
optics, theory of water waves ^and in other areas. 

To fix ideas we assume the potentials to be of the form Vh{x) := V{hx) 
with V smooth and decaying at oo. Thus for h — 0, Equation (QJ becomes the 
standard generalized nonlinear Schrodinger equation (gNLS) 

^^^-Ai; + ^li;-fm')^, (2) 

*This paper is part of the first author's Ph.D thesis. 

t Supported by NSERC under Grant NA7901 and NSF under Grant DMS-0400526. 
^In this case Equation (0 is called the Gross-Pitaevskii equation. 

^In these two areas Vti arises if one takes into account impurities and/or variations in 
geometry of the medium and is, in general, time-dependent. 



where = V{0). For a certain class of nonlinearities, /(IV'P) (see Section (S)), 
there is an interval Xq C K such that for any X G Iq Equation ( |2) has solutions 
of the form e'^^' (j)^ {x) where 0^ G ?i:2(M") and (j)^ > 0. Such solutions (in 
general without the restriction 0q > 0) are called the solitary waves or solitons 
or, to emphasize the property 0q > 0, the ground states. For brevity we will use 
the term soliton applying it also to the function (j)^ without the phase factor 

gi(A-/j)t^ 

Equation ([^J is translationally and gauge invariant. Hence if e*^'^~^-'*0o (^) 
is a solution for Equation ([5J), then so is 

g^(A-M)te«"(^^(2; + a), for any a e K", and a e [0, 27r). 

This situation changes dramatically when the potential Vh is turned on. In 
general, as was shown in jFW Ohl out of the (n + 2)-parameter family 
e*^'^~^-"e'"0Q(a; + a) only a discrete set of two-parameter families of solutions to 
Equation (P) bifurcate: e*'^*e'"0'^(a;), a € [0, 27r) and A G I for some I C Ig, 
with (f)^ = (p^ € 7^2 (K") and 4>^ > Q. Each such family centers near a different 
critical point of the potential Vh{x). It was shown in |Oh2| that the solutions 
corresponding to minima of Vh {x) are orbitally (Lyapunov) stable and to max- 
ima, orbitally unstable. We call the solitary wave solutions described above 
which correspond to the minima of Vh{x) trapped solitons or just solitons of 
Equation ( ^ omitting the last qualifier if it is clear which equation we are 
dealing with. 

Results. In this note we describe results of |GS1I IGS2| that the trapped 
solitons of Equation ( ^ are asymptotically stable. The latter property means 
that if an initial condition of ( ^ is sufficiently close to a trapped soliton then 
the solution converges (relaxes), 

V'(a;,i) -e*'^(*V^°°(2:) ^ 0, 

in some weighted space to, in general, another trapped soliton of the same 
two-parameter family. We also find effective equations for the soliton center 
and other parameters. In this paper we prove this result under the additional 
assumption that if d > 2 then the potential is spherically symmetric and that 
the initial condition symmetric with respect to permutations of the coordinates. 
In this case the soliton relaxes to the ground state along the radial direction. 
This limits the number of technical difficulties we have to deal with. We expect 
that our techniques extend to the general case when the soliton spirals toward 
its equilibrium. 

In fact, |GS1I IGS2j prove a result more general than asymptotic stability 
of trapped solitons. Namely, we show that if an initial condition is close (in 
the weighted norm := ||(1 -I- for sufficiently large u) to the 

soliton e^'">(j)^°, with 70 G M and Aq G X (X as above), then the solution, ip{t), 
of Equation ( ^ can be written as 

V'(x, t) = e*T(*) (e*P(*)-"(/)^(*) {x - ait)) + R{x, t)) , (3) 
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where 0, X{t) Aoo for some Aoo as f — + cxj and the sohton center 

a{t) and momentum p{t) evolve according to an effective equations of motion 
close to Newton's equation in the potential h'^V{a). 

We observe that ( ^ is a Hamiltonian system with conserved energy (see 
Section ^ and, though orbital (Lyapunov) stability is expected, the asymptotic 
stability is a subtle matter. To have asymptotic stability the system should be 
able to dispose of excess of its energy, in our case, by radiating it to infinity. 
The infinite dimensionality of a Hamiltonian system in question plays a crucial 
role here. This phenomenon as well as a general class of classical and quantum 
relaxation problems was pointed out by J. Frohlich and T. Spencer [Private 
Communication] . 

We also mention that because of slow time-decay of the linearized propaga- 
tor, the low dimensions d — 1,2 arc harder to handle than the higher dimensions, 
d>2. 

Previous results. We refer to jGSl| for a detailed review of the related 
literature. Here we only mention resuhs of 'Uu "BPl BP2 BISul IgWTl [5W2l 
I5W31 ITTTl ITY2I It73] which deal with a similar problem. Like our work, 
|Wlll?W2l[nW8llT711lT72llT73l study the ground state of the NLS with a 
potential. However, these papers deal with the near-linear regime in which the 
nonlinear ground state is a bifurcation of the ground state for the corresponding 
Schrodinger operator —A + V{x). The present paper covers highly nonlinear 
regime in which the ground state is produced by the nonlinearity (our analysis 
simplifies considerably in the near-linear case). Now, papers |Cu[ IBPll IBP2L 
IBuSuj consider the NLS without a potential so the corresponding solitons, which 
were described above, are affected only by a perturbation of the initial conditions 
which disperses with time leaving them free. While in our case they, in addition, 
are under the influence of the potential and they relax to an equilibrium state 
near a local minimum of the potential. 

Open problems. We formulate some open problems: 

(1) Extend the results of this present paper to more general initial conditions 
and to more general, probably time-dependent, potentials. 

(2) Link the results of this paper with the results of jFGJS| on the long 
time dynamics of solitons. 

A natural place to start here is spherically symmetric potentials but general 
initial conditions. Note that for certain time-dependent potentials the solitons 
will never settle in the ground state. 

Notation. As customary we often denote derivatives by subindices as in 
^\ ~ Wx't'^ ^'^'^ 't'^ ~ 4>^{^)- However, the subindex h signifies always the de- 
pendence on the parameter h and not the derivatives in h. The Sobolev and 
spaces are denoted by Ti,^ and respectively. 

Acknovirledgment. We are grateful to J. CoUiander, S. Cuccagna, S. Dejak, 
J. Frohhch, Z. Hu, W. Schlag, A. Soffer, G. Zhang, V. Vougalter and, especially, 
V.S. Buslacv for fruitful discussions. This paper is part of the first author's Ph.D 
thesis requirement. 
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2 Hamiltonian Structure and GWP 



Equation is a Hamiltonian system on Sobolev space 7i^(R, C) viewed as a 
real space n^{m, M) ® (R, M) with the inner product {ip, (j)) = Re /jj "00 and 
with the simpletic form uj{ip,(f)) — Im Jj^tp4'- The Hamiltonian functional is: 

HW ■■= J[li\^.\' + V,m')-F{\^pr)], 

where F{u) 

Equation ( ^| has the time-translational and gauge symmetries which imply 
the following conservation laws: for any i > 0, we have 

(CE) conservation of energy: H{'4){t)) = H{ip{0)); 

(CP) conservation of the number of particles: N{tlj{t)) = N{ip (0)), where N{ip) := 

To address the global well-posedness of ( ^ we need the following condition on 
the nonlinearity /. Below, s+ = s if s > and = if s < 0. 

(fA) The nonlinearity / satisfies the estimate |/ (^)| < c(l + |^|"^^) for some 
" e [0, and 1/(01 < c(l + lei'^) for some /3 e [0, |). 

The following result can be found in |Cazj . 
Theorem Assume that the nonlinearity f satisfies the condition (fA), and that 
the potential V is bounded. Then Equation is globally well posed in Ti} , i.e. 
the Cauchy problem for Equation with a datum 0(0) G Tl^ has a unique 
solution ip{t) in the space Ti.^ and this solution depends continuously on "0(0). 
Moreover ip{t) satisfies the conservation laws (CE) and (CP). 



3 Existence and Orbital Stability of Solitons 

In this section we review the question of existence of the solitons (ground states) 
for Equation ( ^ . Assume the nonlinearity / : M ^ R is smooth and satisfies 

(fB) There exists an interval Xq G M+ s.t. for any A G lo, —oo < lims^+oo — — < 
and i f{s)ds > A for some constant ^, for d > 2; and 

C/(0,A) :=-A02+ / /(Ode 

has a smallest positive root 0o(-^) such that U^{(j)o{^), A) 7^ 0, for d = 1. 

It is shown in |BLI IStr| that under Condition (fB) there exists a spherical 
symmetric positive solution 0^ to the equation 

-A0^ + A0^-/((0^)2)0^=O. (4) 
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Remark 1. Existence of soliton functions (j)^ for d — 2 is proved in JA'^jf under 
different conditions on f . 

When the potential V is present, then some of the sohtons above bifurcate 
into sohtons for Equation (^. Namely, let, in addition, / satisfy the condition 
1/ (01 — '^(1 + 1"^!^); foi' some p < oo, and V satisfy the condition 

(VA) V is smooth and is a non-degenerate local minimum of V. 

Then, similarly as in )FW[ IOhl| one can show that if h is sufficiently small, 
then for any A G Iqv, where 

Jov {A|A > - M{V{x)}} n {A|A + ^(0) e Xo}, 

there exists a unique soliton (p^ = 0^ (i. e. 0'^ e 7^2 (K) and > 0) satisfying 
the equation 

-A^^ + (A + - fii<j)^)^)^^ = 

and the estimate 0'^ = cj)'^^^'"^'' + 0{h^/^) where 0q is the soliton of Equation 
(II- 

Let (5(A) :— H^'^Hi- It is shown in |GSSlj that the soliton (f)'^ is a minimizer 
of the energy functional H{ip) for a fixed number of particles N{ip) = constant 
if and only if 6 (A) > 0. Moreover, it is shown in |We2l IGSSII that under 
the latter condition the solitary wave ^^e*'*'* is orbitally stable. Under more 
restrictive conditions (see IGSSlp on / one can show that the open set 

X := {A e lov ■■ S'{X) > 0} (5) 

is non-empty. Instead of formulating these conditions we assume in what follows 
that the open set T is non-empty and X E T. 

Using the equation for cf)^ one can show that if the potential V is re- 
dially symmetric then there exist constants c, 6 > such that |(/)'^(x)| < 
(.g-s\x\ g^j-^j < ce^'^l^l, and similarly for the derivatives of (j)^ and -^4>^- 

4 Linearized Equation and Resonances 

We rewrite Equation (^ as ^ = G{ip) where the nonlinear map G{ip) is defined 
by G{ip) = — i(— A + A + Vii)tlj + if{\ip\'^)ip. Then the linearization of Equation 
(nj can be written as ^ = dG{(f)^)x where dG{(j)'^) is the Frechet derivative of 
G{tp) at (f). It is computed to be 

dG{<j>^)x = -i{-A + X + Vh)x + if({<l^^?)x + 2if\{cp^f){<l,^fRex. (6) 

This is a real linear but not complex linear operator. To convert it to a linear 
operator we pass from complex functions to real vector-functions x ^ — ^ X = 
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where xi = R^X ^J^d X2 = Ittix- Then dG{(p^)x < — > L(yX)x where 



V X2 

the operator L{X) is given by 



with L_(A) -A + 14 + A-/((0^)2), and L+(A) _A + F;, + A - - 
2/'((0^)2)(0^)2. The operator L(A) is extended to the complex space '^^(R, C)© 
7i^(]R, C). If the potential Vh in Equation decays at oo, then by a general 
result 

aess{L{X)) ~ {—ioo, —iX] n [iX, zoo). 
The eigenfunctions of L{X) are described in the following theorem (cf |GS1| . 

Theorem 4.1. Let V satisfy Condition (VA ) and \h\ be sufficiently small. Then 
the operator L{X) has at least 2d + 2 eigenvectors and associated eigenvectors 
with eigenvalues near zero: two-dimensional space with the eigenvalue and a 
2d- dimensional space with non-zero imaginary eigenvalues ±iej(A), 

e^(A) := h^/2e] + o{h), 



where Cj are eigenvalues of the Hessian matrix ofV at value x = 0, V (0). The 
corresponding eigenfunctions ^ ^j,^. ^ '^'^^ related by complex conjugation and 
satisfy 

— \/29a;j,0Q + o(/i) and r]j = —hy/ejXj<pQ + o{h), 
and and rjj are real. 

Remark 2. The zero eigenvector { \\] o,nd the associated zero eigenvector 
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are related to the gauge symmetry ipix, t) — > e'"'0(a;, t) of the original 
equation and the 2d eigenvectors ^ . ^ with 0(h) eigenvalues originate 
from the zero eigenvectors ( ° 1 , fc = 1, 2, • • •, d, and the associated zero 

eigenvectors ^ ^ ^ , k ~ 1,2,- • -jd, of the = equation due to the 

translational symmetry and to the boost transformation ip{x,t) e^^'''ip{x,t) 
(coming from the Galilean symmetry), respectively. 

For d > 2 we will be interested in permutationally symmetric functions, 
g e £^(M''), characterized as 

g{x) — g{crx) for any a ^ Sd 

with Sd being the group of permutation of d indices and a{xi,X2,- ■ -^Xd) '■= 

(2^(7(1), 2:cr(2), • • -iXa^d))- 
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Remark 3. For any Junction of the form e^^'^(j){\x — a\) with a \\ p, there exists 
a rotation r such that the function e^P"^^(j){\TX — a\) = e'"^ P'^(j){\x — T~^a\) 
is permutationally symmetric. Such families describe wave packets with the 
momenta directed toward or away from the origin. 

If for (i > 2 the potential V{x) is spherically symmetric, then V (0) = 
4Ay(0) • Idnxni and therefore the eigenvalues Cj of V (0) are all equal to 
|Ay(0). Thus we have 

Corollary 4.2. Let d>2 and V satisfy Condition (VA ) and let V he spherically 
symmetric. Then L{X) restricted to permutational symmetric functions has 4 

eigenvectors or associated eigenvectors near zero: two-dimensional space with 
eigenvalue 0; and two-dimensional space with the non-zero imaginary eigenval- 
ues ±ie(A), where 

and with the eigenf unctions ( ^.^ ], where ^ and t] are real, and permuta- 



±iri{X) 

tion symmetric functions satisfying 



V2 ^ -^(/.^ + 0{h) and r?(A) = ~hJ^AV{0) ^ + ©(/i^/^). 

dXn \ d ^ 

n=l n=l 



The eigenvectors ( > s, ) are symmetric combinations of the eigenvec- 



tors described in Theorem 2. 

Besides eigenvalues, the operator L(A) may have resonances at the tips, ±iA, 
of its essential spectrum (those tips are called thresholds). Recall the notation 
a+ := a if a > and = of a < 0. 

Definition 4.3. Let d ^ 2. A function h is called a resonance function of L{X) 
at jjL = ±iA if h ^ L^, \h{x)\ < c(x)~^'^~^^+ and h is and solves the equation 

(L(A) - n)h = 0. 

Note that this definition implies that for d > 2 the resonance function h 
solves the equation {l+K{X))h = where K{X) is a family of compact operators 

given by if (A) := {Lo{X)- n+OyWugW- HereLo(A) ■■=(^\_x 0^^^) 
and 

'''"'^>-\-VK+f{{0'r)+v'{{^?){^? )' ^ ^ 

In this paper we make the following assumptions on the point spectrum and 
resonances of the operator L(A) : 
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(SA) L{X) has only 4 standard and associated eigenvectors in the permutation 
symmetric subspace. 

(SB) L{X) has no resonances at ±iA. 

The discussion and results concerning these conditions, given in 'GSl', sug- 
gested strongly that Condition (SA) is satisfied for a large class of nonlinearities 
and potentials and Condition (SB) is satisfied generically. Elsewhere we show 
this using earlier results of jCPI ICPVj . We also assume the following condition 

(FOR) Let N be the smallest positive integer such that e(A)(iV + 1) > A, VA e /. 
Then ReYjs; < where y„, n = 1,2, ■ ■ ■ , are the functions of V and A, 
defined in Equations f [T7|l below (see also (El- 

We expect that Condition (PGR) holds generically. Theorem 15.21 below 
shows that ReYn ^ Q ii n < N. 

We expect the following is true: 
(a) if for some Ni{> N), ReYn = for n < Ni, then RcYn, < and (b) for 
generic potentials/nonlinearities there exists an Ni{> N) such that ReY^-^ ^ 0. 
Thus Condition (PGR) could have been generalized by assuming that RcYn^ < 
for some Ni > N such that ReYn = ior n < Ni. We took = A^i in order 
not to complicate the exposition. 

The following form of ReY^ 

ReYN = Im{(Ji{L{\) ^ {N + l)ie(A) - Oy^F, F) < (9) 

for some function F depending on A and V and ai :— ^ , is proved 

in iTYn lTY2l lTY8l [SW4j for A^ = 1, and in fej for A^ 2, 3. We 

conjecture that this formula holds for any A^. 

Condition (PGR) is related to the Permi Golden Rule condition which ap- 
pears whenever time-(quasi)periodic, spatially localized solutions become cou- 
pled to radiation. In the standard case it says that this coupling is effective 
in the second order (A^ = 1) of the perturbation theory and therefore it leads 
to instability of such solutions. In our case these time-periodic solutions are 
stationary solutions 

c,(fJ."«' + c.(l,Je-.«' 

of the linearized equation ^ — L{X)x and the coupling is realized through 
the nonlinearity. Since the radiation in our case is "massive"— the essential 
spectrum of L{X) has the gap {—iX,iX), A > 0, — the coupling occurs only in 
the A^— th order of perturbation theory where A^ is given in Condition (PGR). 

The rigorous form of the Permi Golden Rule for the linear Schrodinger equa- 
tions was introduced in |BS| . Por nonlinear waves and Schrodinger equations 
the Permi Golden Rule and the corresponding condition were introduced in |S] 
and, in the present context, in |SW4l IBuSul lBP2l ITYTI ITY21 ITYS] . 
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5 Main Results 

In this section we state the main theorem of this paper. For technical reason 
we impose the fohowing conditions on / and V 

(fC) the nonUnearity / is a smooth function satisfying / (0) = / (0) = if 
d>3; and /(''^(O) = for fc = 2, 3 • • • 37V + 1 if d = 1, where /('^^ is the 
fc— th derivative of /, and N is the same as in Condition (FGR), 

(VB) V decays exponentiaUy fast at oo. 

Theorem 5.1. Let Conditions (JA)-(fC), (VA), (VB), (SA), (SB) and (FGR) 
be satisfied and let, for d > 3, the potential V be spherically symmetric. Let an 
initial condition V'o be permutation symmetric if d > 5 and A G X. There exists 
c, eo > such that, if 

(10) 

with \{zi,Z2)\ < eq and n — 1,2 being real, some large constant v > and 
with k = [^] + 2 if d > 3, and k ~ 1 if d ^ 1, then there exist differentiable 
functions 7, zi, Z2 : M"*" ^ R, A : IR+ I and R : TL^ such that the 

solution, ijj{t), to Equation ('QJ) is of the form 

^{t) = e'^o Ks)ds+^^{t)^^X{t) ^ ^ ^ 

with the following estimates: 

(A) II (1 + a;^)^''i?(t)||2 < c(l + " where V and N are the same as that in 
( \1U\) and (FGR) respectively, 

2 

(B) Y,\z,{t)\<c{l + t)-^. 
i=i 

Remark 4. Recall from Remark that the class of permutationally symmet- 
ric data includes wave packets with initial momenta directed toward or in the 
opposite direction of the origin. 

Theorem 5.2. Under the conditions of Theorem 3 we have 

(A) there exists a constant Aqo G T such that hm \{t) = Aqo- 

t — >co 

(B) Let z := zi — iz2. Then there exists a change of variables (3 = z + 0{\z\^) 
such that 

N 

(5 = te{X)(3 + J2 Yn{X)P''+'P" + 0(|/3|2~+2) (^2) 

n=l 

with Yn being purely imaginary if n < N and, by Condition (FGR) ReY^ < 0. 
Moreover, for N — 1,2,3, ReY^ is given by Equation (\^. 
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Remark 5. Using that e(A) = h^^^^ + o{h) one can rewrite Equations f \llp 

and ( liijj) in the form with a(t) and pit) satisfying the equations = p 
andp^ -h^VV{a) modulo 0(|ap + 

6 Idea of the Proof: Expansions 

We follow |miGS2| . We decompose of the solution ip{t) to Equation into a 
solitonic component and a simplectically orthogonal fluctuation as (cf iFGUS p 

^{t) = e'-lo^^'')'''+'^('\c^^ + zi{t)^ + iZ2{t)rj + R{t)), (13) 

where A, 7, zi, Z2 are real, differentiable functions of t and the function R{t), 
called the fluctuation^ satisfies the orthogonality conditions 

Im{R, i(j)^) = Im{R, -^<j)^) = Im{R, irj) Im{R, = 0. (14) 

We plug Equation ( I13|l into Equation ( ^| to obtain equations for the pa- 
rameters zi{t), Z2{t), X{t) and 7(t) and the fluctuation R{t). As was already 
discussed above the linearized operator, G{4>^) (see (EJ), in the equation for 
R{t) is only real-linear and therefore we pass from the unknown R to the un- 

known i? := ( \ ^ R. Under this correspondence the multiplication by 

goes over to the symplectic matrix (^^ ^ ) ' ^ i^^R- Unlike 

with the equation for R, in the equations for zi and Z2 it is more convenient to 
go from the real, symplectic structure given by J to the complex structure i^^ 

by passing from ^ ^ to z := zi — iz2. 

A key point is to look for the fluctuation R{t) in the form ( [nilGS2| ) 

i?= i?™,„(A)z™z" + i?jv (15) 

2<m+ri<Af 

with the remainder Rn of the order 0(|z|^+^). This leads to the following 
equations on the coefficients Rmn- 

[L(A) - ie{\){m - n)]4„„(A) = -Pc/m,„(A), 

where the functions fm.nW depend on Rm'.n'W with m' + n' < m + n. Recall 
that if |r7i — 7i| < N, then i£{X){m — n) ^ (t(L(A)) and therefore the operators 

L{X) ^ ie{X){m - n) : P,,C'^ ^ PcC"^ (16) 

are invertible. Hence the above equations have unique solutions. 

We plug the expansion ( 1151) the differential equations for the parameters 
A, 7 and z to obtain expansions for A, 7 and z in terms of z and z. 
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In the second step we transform z to a parameter y which satisfies a simpler 
different equation. We construct a polynomial P{z, z) with real coefficients and 
the smallest degree > 2, such that the new parameter y:= z-\- P{z, z) satisfies 
the equation 

y = ie{\)y+ ^ Y^^nWy^r 

2<m+n<2JV+l (17) 
^Remainder. 

where the coefficients Ym,n{^) are purely imaginary, and Ym,n{^) = if m 7^ 
n + 1, and the term Remainder admits the estimate 

\Remainder\ < cr'-^ or + \y{t)\^\\{x)-''RN\\2 + \\{x)-''Rn\\1 

Expressing the variables z and z as power series in y and y and plugging 
the result into the z— expansions for 7, A and R mentioned above we obtain the 
y— expansions for these quantities. 

As was mentioned above, the decrease of z (or y) and therefore the relaxation 
of the soliton to its equilibrium occurs due to the radiation of the excess of 
the energy to infinity. The latter is possible only if the periodic solutions to 
the linearized equation are coupled to its continuous spectrum. To detect this 
coupling we must obtain an expansion of R in the parameters y and y up to the 
order 2^^. This is done in the third, most involved step. As in the first step 
we determine the coefficients i?rn,n(A) of the y— expansion of R by solving the 
equations 

[L(A) - ie{X){m - n)]i?„,„(A) = -Pc/m,n(A) 

for certain functions fm,nW (see below). Recall that the number N is defined 
by the properties 

ie(X)(m-n) ^ "^^^^^^ if 1™ - n| < AT, 
ie{A){m n) ^ ^^^x)) if \m-n\> N. 

Thus for N < m + n < 2N the parameter ie(X)(m — n) might be in the spectrum 
of L{X). To deal with this case we sort out the pairs {m, n) into "non-resonant 
pairs" satisfying \m — n\ < N and "resonant pairs" satisfying |m — n| > N. For 
"non-resonant" pairs the operators 

L(A) - ie{X){m - n) : PcC^ P^C^ 

are invertible and for resonant pairs they are not (one has to change spaces in 
the latter case). 

In the first two steps we expanded in z and z (and in y and y) until m+n < N 
and consequently all the pairs, (m, n), involved were non-resonant ones. Now our 
expansion involves pairs (m, n) with m + n > N, which include resonant pairs. 
A key point which we show is that for the subsets of pairs (m, n), m + n > N, 
determined by the inequality 

m,n< N, 
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the terms /m,n(A) involve only "non-resonant" pairs (here we use that the pa- 
rameter y satisfies i Wj\ with F„i,n(A) — Q iov m ^ n + 1) and consequently we 
are able to solve for the coefficients Rm,nW also in this case. 

Finally, estimates of the remainders as in Equations f ll5|) are done by rewrit- 
ing the differential equations for them in an integral form (using the Duhamel 
principle) and using estimates of linear propagators derived in [BP 11 IBP2L 

iE^iH^[nH;^[n!?Ti[nn2] . 
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